Speeding up the spatial adiabatic passage of matter waves in optical microtraps by 

optimal control 
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We numerically investigate the performance of atomic transport in optical microtraps via the so 
called spatial adiabatic passage technique. Our analysis is carried out by means of optimal control 
methods, which enable us to determine suitable transport control pulses. We investigate the ultimate 
limits of the optimal control in speeding up the transport process in a triple well configuration for 
both a single atomic wave packet and a Bose-Einstein condensate within a regime of experimental 
parameters achievable with current optical technology. 
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I. INTRODUCTION 



The coherent control of matter waves has become a 
very relevant research topic with significant technologi- 
cal applications such as in atom lasers |lH3i] and quan- 
tum information processing (QIP) [1], to only cite a few. 
Indeed, very recently, experimental demonstrations with 
ultracold atoms in two-dimensional (2D) optical lattices, 
of single quantum bit (qubit) rotations [5|, single site 
addressability 0], and single atom detection @, 
with very high fidelities, have been reported. One of 
the current most important goals of QIP is to go be- 
yond the manipulation of a handful of qubits since the 
main challenge is to build scalable quantum hardware 
where several thousands of qubits are coherently manip- 
ulated within the relaxation times of the system jlO| . 
It has been also understood that the practical realiza- 
tion of QIP requires to devise new architectures. Indeed, 
several schemes for two-qubit quantum gates implemen- 
tations based on atomic systems proposed in the past 
can be, at least in principle, realized in experiments, but 
they may present several shortcomings when building a 
scalable quantum computing hardware (see Ref. [T]| for 
a review on QIP with neutral particles). Moreover, in 
order to name a QIP system "scalable" , it is also impor- 
tant that the resources required to control the quantum 
system, typically classical devices (e.g., laser fields, re- 
frigerators) , are scalable as well [l^l ■ 

Paradigmatic examples of QIP are the Cirac-Zoller [l^ 
and the M0lmer-S0rensen |14l | ion quantum computers, 
which have been proven to be powerful schemes to ex- 
perimentally realize small quantum algorithms [TBI - flTj 
or to engineer quite exotic entangled states (up to 14 
qubits), like Werner [31 oi Greenberger-Horne-Zeilinger 
states [l^, [l^l- These schemes, however, present rather 
difficult technical problems when hundreds or even thou- 
sands of ions participate in the collective motion (e.g, 
decoherence of motional modes, sensitivity to electric 
noise [2l|). Thus, currently, a big effort is made in the 
design and practical realization of new schemes that are 
actually scalable. For instance, in Ref. [12] it has been 



proposed an ion quantum processor architecture, where 
some areas of the chip processor are used only to store 
the information, and others to manipulate it. Such a 
design requires to transport an ion from one location 
to another one in the chip preserving its quantum me- 
chanical coherence. A similar problem is encountered for 
QIP implementations with neutral atoms either in opti- 
cal lattices [23l - [25| or in microwave atom chips [1^ [l^] . 
Atoms trapped in optical lattices can be efficiently pre- 
pared via the superfluid Mott insulator quanturnphase 
transition [2^, [2^ , and single sites addressed [^, [J, but 
the realization of quantum gates between qubits located 
in far away lattice sites can be a serious problem for a 
scalable neutral-atom-based quantum processor. Solu- 
tions to this issue can be afforded, for instance, by aux- 
iliary atoms that can be efficiently transported in state- 
independent periodic external traps (soj or by using op- 
tical tweezers [3l| . 

A major underlying concept of these QIP paradigms 
is the coherent transport of ions or atoms in such a way 
that they mediate the operation of quantum gates be- 
tween spatially distant qubits. The needed transport 
time, however, has been estimated to be about 95% of 
the time used for carrying out the whole quantum com- 
putation [32j . It is therefore imperative to reduce the 
time needed to transport an atom or ion from the quan- 
tum memory to the processing units and, therefore, to 
engineer robust control transport pulses. In this respect, 
optimal control theory is a prominent candidate for a 
drastic improvement of the design of accurate QIP pro- 
tocols, and, recently, several theoretical investigations 
on the optimal transport of both a sin gle atom and an 
atomic ensemble have been undertaken [33-|3g. Besides 
this, very recently, control pulses numerically obtained 
by using iterative optimization algorithms have been ex- 
perimentally applied, with great success, in order to effi- 
ciently transfer a one-dimensional (ID) degenerate Bose 
gas from the transverse ground to the lower excited state 
of a waveguide potential [s^- This result shows the po- 
tential offered by optimal control methods to engineer 
current experiments of ultra-cold atoms. 

In this work we investigate the ultimate limits of the 
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transport of neutral atoms in optical microtrap arrays by 
means of numerical optimization methods. Specifically, 
we are interested in the spatial adiabatic passage (SAP) 
protocol [13 , the matter wave analogue of the stimulated 
Raman adiabatic passage (STIRAP) technique used in 
quantum optics to transfer the population between two 
atomic internal levels [s^- The SAP technique consists 
in adiabatically following an energy eigcnstate of the sys- 
tem, the so-called spatial dark state, that only involves 
the vibrational ground states of the two extreme wells 
of a triple- well potential (see Fig. [T]). The spatial dark 
state presents at all times a node in the central region 
such that the middle-well population is almost negligi- 
ble throughout the transport process. Therefore, the 
SAP protocol enables to transport an atom from a lat- 
tice site to the next-nearest-neighbor without populating 
the nearest-neighbor site (the middle well in Fig. [T]). If 
one atom is present in the middle trap, the SAP tech- 
nique can be used to implement a single atom diode or 
a transistor (40| . In addition, this transport technique 
may allow to reduce the complexity of several quantum 
computing architectures [tI. [26l. l4l| - |43| . In fact, compared 
to the tunneling induced oscillation between two adjacent 
traps, such technique is more robust against variations of 
the system parameters and requires less precise control 
of the distance and timing. Three-level atom optics tech- 
niques [13 , such as the SAP protocol, allow also to create 
superpositions (spatial dark states) of matter waves be- 
tween two separated sites of the lattice, useful for appli- 
cations in atomic interfcromctry, or to inhibit the tunnel- 
ing among lattice sites, therefore allowing to create con- 
ditional phase shifts for quantum logic, or to transport 
an empty site [s^- We also mention, that, recently, the 
implementation of the SAP protocol for radio-frequency 
traps [43 - l46j within the three mode approximation has 
been investigated [i^l, and that such technique could be 
even employed for the ion transport in segmented micro- 
traps |48l. l49j. 
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Figure 1. (Color online). Upper panel: Initial (black thick 
line) and final (blue thin line) potential configuration for step 
1 of the transport process with trap separation 6.5 /im, beam 
waist 2w = 1.3 /^m, potential depth Vb = fcs x 86 uK, and 
trap frequency ojx — 4Vo jma'^ — 2ttUx ~ 27r x 711 Hz. 
The initial ground state wave function i/'o(a;) = 4i^^l{x) (red) 
of the left well is also displayed. Lower panel: Initial and 
final potential configuration (same as the blue thin line of 
the upper panel) for the SAP process. The goal wave func- 
tion i'gix) = (l>'i^l{x) (bright grey) of step 1 is shown. This 
state also corresponds to the initial wave function of the SAP 
process, whose goal wave function is ^g{x) = <j>^^\x) (dark 
green), superimposed on the right well of the lower panel. 
The horizontal black arrows indicate the transport direction. 



II. OPTICAL DIPOLE MICROTRAPS 



Beside this, we will investigate the transport of a Bose- 
Einstein condensate (BEC) in relation to recent experi- 
ments with optical dipole traps [sol - fs^ . We note that a 
similar study has been carried out in Ref . [33| , where the 
optimal transport of a BEC in magnetic microtraps, like 
the ones produced with atom chips [1], has been investi- 
gated, and that, very recently, the optimal control pulses 
for harmonically trapped BECs have been analytically 
determined [1^. We underscore that, while the goal of 
those investigations was to transfer a BEC between spa- 
tially separated locations, here, in addition to this goal, 
we aim at minimizing the occupancy of the middle well 
in a triple-well configuration, as showed in Fig. [l] by 
following as much as possible the spatial dark state of 
the trapping potential. This additional constraint is the 
main signature of the SAP protocol. 



We consider a (transverse) potential, where either a 
single atom or BEC is trapped, given by the following 
analytical expression 
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where Vq represents the depth of the three Gaussian 
dipole traps, with 2w being the laser beam waist, d^iit) 
represents the distance between the central trap and the 
left trap, di{t) is the distance between the central trap 
and the right trap, while the central trap remains at 
rfo(t) = Vt e [0,T]. Here T is the time needed to 
transport the system of interest (i.e., an atom or a BEC) 
initially prepared in the ground state of the trap on the 
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left (centered in — I^qI) to the ground state of the trap 
on the right (centered in \xq\). As outhned above, we 
shall first consider a ID scenario, but later we shall also 
study the influence of the dimensionality on the transport 
performance. We note that optical dipole traps, as the 
ones of Ref. 
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can be designed to form a 2D lattice. 
The X — y plane, that defines the lattice, has a weaker 
confinement than the (vertical) axial direction, thereby 
defining a "pancake" geometry. It is in the (transverse) 
x — y plane that the transport occurs and because of this 
we refer to the potential ([T]) along x as transverse (in our 
study the trap frequencies in the x and y directions will 
be assumed to be equal). 

In Fig.[T]thc potential is illustrated, where experimen- 
tally realistic parameters (i.e., potential depth, trap sep- 
aration, and beam waist) have been considered for ®^Rb 
atoms. Due to the large (initial) separation between the 
trap minima of the lattice [x'q = 6.5 /xm, see Fig. [1] (top)], 
no tunneling is expected to happen until the traps are 
closer. The SAP transport is split and optimized in three 
different stages. Firstly, in a time Ti, the initial atomic 
state (the ground state of the left well) is moved from 
x'q ~ —6.5 fim to a;o = —2.5 fiui, i.e., from the left well of 
the initial trap configuration (thick black line) given in 
Fig. [1] (top) to the left well of the potential of the upper 
panel of Fig. [1] (blue thin line). Secondly, the atomic sys- 
tem is brought, in a time T2, from the left well centered at 
xq ~ —2.5 fim to the right one centered at Xq = 2.5 /im 
(see blue line of the lower panel of Fig. [T|). The third 
step consists in bringing the system from the right well 
centered at xq = 2.5 fim to the right well centered at 
x'q = —6.5 /im of the potential displayed in Fig. [T] (top), 
which is equivalent to reversing the first step of the pro- 
cess. Hence, the total transport time is T = 2Ti + T2. It 
is in the second step that the SAP process takes place. 
Since at the beginning the atoms are quite far apart, 
the direct application of the SAP technique would sim- 
ply slow down the whole transport process, because ini- 
tially no tunneling would take place. Instead, with the 
above outlined procedure, the first step can be sped up 
as much as possible until the so-called quantum speed 
limit is reached, that is, the minimum time required to 
evolve a quantum system from an initial state to an other 
orthogonal state |53l - [5a |. 

Finally, we note that a crucial condition for the re- 
alization of SAP is that the initial and final states in- 
volved during the transport process should be in reso- 
nance. This requirement can be fulfilled by fixing at all 
times the minima of the potential at the same energy 
level, through the control of the time-dependent param- 
eters V-i{t), vi{t) [vo{t) = I] (sec also Fig. [2|), as well 
as by fixing the maxima of the triple well configuration 
to the same energy level, through the control of w, that 
is, the beam waist (see the appendix for the analytical 
expressions of v±i). The control of the latter, however, 
would require beam waist values below the actual experi- 
mentally achievable limit (sij , and therefore it will not be 
considered in our study. Hence, in our analysis we have 
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Figure 2. (Color online). Potentials at time t = T2/2: 
the blue (solid) line shows the potential when v±i are time- 
dependent, whereas the magenta (dashed) line when v±i are 
time-independent. The potential minima have been lifted by 
about 0.33/i;/a; and OAihVx, respectively. The symmetry of 
the potentials is due to the fact that, at that time, the two 
outer traps are equidistant from the centre a; = 0, that is, a 
distance almost equal to the minimal allowed trap separation 
Sxq = 1.43 fim. The rest of trap parameters are as in Fig. [T] 



fixed w to the minimum experimentally achievable value 
(i.e., 0.65 /im), which enables us to prepare the atoms, 
before the transport, in a lattice configuration with min- 
imum periodicity. 



III. SINGLE ATOM TRANSPORT 

In this section we analyze the transport of a single 
trapped atom in the absence of a thermal or quantum 
bath. The atomic state obeys the Schrodinger equation 
of motion. In order to speed up the transport process we 
rely on numerical optimization techniques. For the prob- 
lem considered in this paper, we employed a recently in- 
troduced optimization method, named the chopped ran- 
dom basis (CRAB) algorithm [HBl- Such a method has 
been shown to be a powerful tool in order to optimize the 
closed dynamics of many-body quantum systems (STj and 
the dynamics of light harvesting (5^. Moreover, since the 
implementation of the CRAB algorithm does not rely on 
the equation of motion that governs the system dynam- 
ics, there is no need of algorithmic modifications when 
nonlinear dynamics is regarded (e.g., the dynamics of a 
BEC), unlike for the monotonically convergent Krotov 



algorithm 159, ,61 
ing algorithm [62 



or the gradient ascent pulse enginecr- 
. Even though CRAB does not provide 
monotonic convergence it allows to directly restrict and 
select the space of control pulses (e.g., enforcing limited 
bandwidth), since it relies on a multi- variable function 
minimization that can be performed, for example, via 
a direct-search method (e.g., the Nelder-Mead method 
as implemented, for instance, in MATLAB). For more 
details on the procedure of the CRAB algorithm imple- 
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mentation and its computational performance we refer to 
Ref. lEi. 



A. Optimization of step 1 of the transport pro 



icess 



We remind that our goal here is to transport an atom 
initially prepared in the ground state of the left well of 
the potential displayed in Fig. [1] (top) , centered in x'q — 
—6.5 iim, to the ground state of the left well of the poten- 
tial shown in Fig. [1] (bottom) blue (solid) line, centered in 
xq ~ — 2.5^m. In this case we shall consider the control 
pulses to be identical, that is, d_i(t) = di{t) = d{t). 

In this section, the objective functional (to be mini- 
mized) for the control problem we are interested in can 
be identified by the so-called overlap infidelity, namely 



da; ■>p'*{x)ip{x,Ti) 



(2) 



Here ip{x, Ti) is the wave function at time t = Ti propa- 
gated from the initial condition iIiq{x) = <j}'}\{x) at time 

t = 0, where 0^^(a;) is the ground state of the left trap 
centered at a; = Xq = —6.5 /im. The wave function il}'g{x) 
is the wave function we aim to achieve in a given time Ti , 

f 2") 

that is, the ground state 4)_{{x) of the left well centered 

aX X ^ Xq = —2.5 /im. The superscript {j) in refers 
to the two first stages of the transport process. When 
j = 1 the state corresponds to the ground state of 
the left well of the potential (thick black line) illustrated 
in the upper panel of Fig.[Tl while for j = 2 it corresponds 
to the ground state of the left well of the potential shown 
in the lower panel of Fig. [T] (blue solid line) . The same 
a) 

applies to <PY , but for the right wells. 

The wells of the upper panel of Fig. [T] (thick and thin 
lines) are sufficiently deep that the trapping potentials 
can be, with good approximation, considered harmonic, 
but with slightly different trap frequencies. Hence, an 
excellent guess control pulse is given by [s^ 



4o(i) = (a;o -xo)\^+ sin 
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(3) 



For a particle in a moving harmonic potential such a con- 
trol pulse is optimal, i.e., it yields 1 = 0, and it is quite 
robust against control pulse distortions [s^l- Neverthe- 
less, as Fig. |3]shows (solid line), for our case with Gaus- 
sian traps, we obtain already a good result for large, but 
not adiabatic (i.e., not in the regime uJxTi ^ 1), trans- 
port times. 

To reduce the infidelity we further optimized the trans- 
port process by means of the CRAB algorithm, which 
works as follows: we start with the initial guess given by 
Eq. and we define the new control pulse as d{t) = 
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Figure 3. (Color online). Overlap infidelity vs. transport 
time: black (solid) line with the control pulse duo{t) defined 
in Eq. ((3|; the other two lines show the infidelity obtained 
with the CRAB optimized control pulse rfopt(i) ~ dho{t)gopt{t) 
for Ng = 8 (dashed line) and Ng — 16 (dot-dashed line). 



duo{t)g{t), where 



1 

lit) = l + TTiylZt^fc sin(c^fct) 
^ ' fe=i 



Bk cos (uJkt)] . (4) 



Here tOk = 27rA:/Ti, iV^ e N is the number of time- 
independent Ak and Bk coefficients, A(t) is a time- 
dependent function enforcing the boundary conditions 
of d{t) at t ~ and t = Ti, namely limt^o.Ti A(t) = 
+00. Basically, the CRAB algorithm seeks for the time- 
independent coefficients Ak, Bk and frequencies uJk thai 
minimize the overlap infidelity Besides, in the nu- 
merical simulations, we set a tolerance (~ 10~^) on the 
determination of either the coefficients or the frequen- 
cies UJk- Such a tolerance is defined as the minimum 
allowed distance between the vertexes of the polytope 
generated within the Nelder-Mead multidimensional non- 
linear minimization procedure (63j . 

As illustrated in Fig. |31 the CRAB algorithm slightly 
improves the result for large Ti times obtained with the 
control pulse defined in Eq. (|3]), but for short times the 
slight difference due to the proximity of the central trap 
in the trap frequencies (of about 0.7%) of the left wells 
[thick and thin lines of Fig. [T] (top)], becomes crucial 
as well as the anharmonicity of the confinement poten- 
tial. We have also investigated the improvement of the 
overlap infidelity due to a higher number of harmonics 
Ng involved in the control pulse. As Fig. [3] shows, a 
significant improvement can be observed only for initial 
large overlap infidelities, whereas for already almost per- 
fect transport processes the reduction is almost insignifi- 
cant. This numerical observation is not surprising, since 
the dynamics with an initial large overlap infidelity re- 
quire more sophisticated control pulses (i.e., higher har- 
monics), in order to properly steer the atomic dynam- 
ics. Such more complex control pulses, however, might 
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Figure 4. (Color online). Optimal control pulses of step 1 of 
the transport process (red dashed line with Ng = 8, blue solid 
line Ng — 16); initial guess given by Eq. ((3| (black thick solid 
line). The transport time is Ti = 3.4 ms. 



Figures. (Color online). Overlap infidelity vs. amplitude of a 
shaking in the trap positions of the outer traps with a;shakG = 
lO^^a;^;, 5x0 = 1.43 /im, and i — \/h/mujx ~ 0.41 pm. The 
transport time is Ti — 3.4 ms and Ng = 16. 



be more difficult to implement experimentally. Beside 
this, as shown in Fig. [U we also note that the overlap 
infidelity drops quite significantly, as a rule of thumb, 
for times larger than 3 ms, which is consistent with the 
fact that the transport time cannot be shorter than the 
inverse of the typical trap frequency, that is, 1/vx ^ 2 
ms. We underscore that the threshold cannot be 
precisely identified with the quantum speed limit, but it 
is in close relation with it. An exact determination of 
the quantum speed limit relics on the time average of 
the instantaneous energy fluctuations, which depend on 
the particular control pulse. This is non trivial compu- 
tational task for time-dependent Hamiltonians. Only in 
simple cases, such as the Landau-Zener model, the min- 
imum time can be efficiently estimated (ssf . However, 
further optimization, that is, larger values of Ng, will 
not overcome this (physical) limit. 

As an example of the optimization. Fig. |4] shows the 
optimized transport control pulses for two different sets 
of coefficients {Ak,Bk} for the transport time Ti — 3.4 
ms, where dho(0 been used as initial guess (black 
thick line). We see that, to achieve very small infideli- 
ties, the control pulse involves more wiggles, especially 
of large amplitude at the intermediate times, where the 
system is more excited. On the other hand, at the end 
of the transport process the system has to be restored 
again in the ground state of the trap, and therefore the 
modulation of the control pulse is more "gentle" . 

We have also investigated the robustness against trap 
position fiuctuations due to possible experimental imper- 
fections. To this aim we used for the time-dependent 
control pulses the following expression 

d^l{t) = dopt{t) ± flshako sin(lLJshakcO) (5) 

where dopt{t) is the optimal control pulse obtained with 
CRAB. With such a choice, for Oshake > the left and 
right wells oscillate in phase, whereas for Oshakc < their 



oscillations are out of phase. However, since in this par- 
ticular step of the transport process only a populated 
well is effectively involved, the only value that matters 
is Iftshakcl (for negative values the behavior of I is basi- 
cally the same) . The result of such analysis is shown in 
Fig. [5] for Wshako/^a; = 10^^ and Ti = 3.4 ms (see also 
Fig. [3]). We see that the optimal control pulse is quite 
robust against fluctuations of the trap position. This re- 
sult is in agreement with the findings for a particle in a 
moving harmonic potential js^ . 

Finally, we investigated the role of spatial dimension- 
ality. Up to now, we performed our analyses in the quasi- 
ID regime. However, we recall that in the experiments 
of Refs. [HO, [EH the potential in the z (axial) direction 
is shallower than in a; or y. We therefore performed nu- 
merical simulations of the 2D Schrodinger equation with 
the trapping potential 

f.^-!^ [x-fcci(t)]^ I 

Here is the beam waist along the z direction. The 
ratio ujz/i-tJx is determined by Wx/wz, namely, the larger 
Wz is, the smaller the (axial) frequency lOz- The new 
(ground) initial and goal states have been obtained by us- 
ing the imaginary time propagation procedure, typically 
used for the determination of the ground state of a BEG. 
As trial functions for the imaginary time propagation we 
took the tensor product of the solutions of the quasi-lD 
regime: for the transverse direction by a numerical exact 
diagonalization of the single particle Hamiltonian, and 
for the axial direction by choosing the Gaussian ground 
state of the harmonic oscillator. 

The result of such analysis, for the optimal control 
pulse obtained for the transport time Ti = 3.4 ms with 
Ng — 16, is shown in Fig. [HI As it is illustrated, 
the smaller uJz is, the better the infidelity. This be- 
havior is reasonable: since Wz S> Wx, we can write 



6 




Figure 6. (Color online). Overlap infidelity vs. the ratio 
ujz/lox- Here X2D is simply the generalization of Eq. ([2]) with 
a two-variable integration. The transport time is Ti = 3.4 ms 
and Ng = 16. 



exp(— 2:^/2zi;^) ~ 1 — /^wl, which implies an almost 
perfect harmonic potential in the axial direction, thereby 
almost separable from the one in the transverse direction. 
The figure shows, however, that the infidelity becomes 
larger than the one obtained in the quasi-lD regime, 
which is almost one order of magnitude smaller (see 
Fig. [3]). We attribute this enhancement to the not com- 
pletely negligible coupling between the axial and trans- 
verse motion. Thus, for an exact infidelity and control 
pulse evaluation, a 2D optimization should be performed. 
Of course, the result relies on the particular chosen uJz/lOx 
ratio, which ultimately is determined by the experimen- 
tal setup. As already pointed out in Ref. [2^ , care has to 
be taken when calculations with realistic parameters in 
the quasi-lD regime are considered. The quantum speed 
limit behavior, however, remains fundamentally the same 
as the one of the quasi- ID regime outlined above, and the 
control pulses obtained in this regime would be excellent 
initial guesses for the 2D optimization. 



B. Optimization of step 2: SAP process 

The transport via SAP is the second step of the opti- 
mization process outlined above, where the initial condi- 

tion at time t = is given by ^o{x) = (f)_l{x), whereas 

(2) 

i^gix) is the ground state wave function (pi (x) of the 
right well centered at x — xq (see Fig.[T]). 

We start the optimization by using the following initial 
guess control pulses: 



(a-o -fao)cos2 



Xo 



5x0 te[0,T2-U] 
te {T2-td,T2] 
(7) 



Xo 



{xo - Sxo) cos^ ( 'rf-fd^ ) + '^^0 



te [0,td] 
te {td,T2]. 



(8) 



We note that di{t) is d^i{t) time inverted. We chose the 
time delay between the two control pulses as <d = 0.17 
where T2 is the transport time used to carry out the SAP 
technique. Such a choice is due to the analogy between 
SAP and STIRAP. Indeed, as shown in Ref. [39|, the time 
over which the two control pulses do overlap, has to fulfill 
the (adiabatic) criterion 



td > 



10 



mintg[o,T]{^(c^(i)/4)}' 



(9) 



where [38 



2^ 



eW^-)'{l + rf[l-erf(d/4)]/4}-l 

^[p2(d/£J2 - 1] 



rf-i(t) 



(10) 

with £x — \/h/moJx- The tunneling "Rabi" frequency de- 
scribes the coupling between the left and the middle wells 
for d = and between the right and the middle wells 
for d = df. We note, however, that Eq. ([TU)) is only valid 
for harmonic trapping potentials. For Gaussian traps the 
actual Rabi frequency must be numerically assessed, but 
for an estimation of the time delay, Eq. ([TOl) can be used. 
As wc will discuss at the end of this section about the 
robustness of SAP against fluctuations of t^, the error in- 
duced by using Eqs. ()9|10p is indeed small, and the value 
of <d used in our analyses is quite reasonable. 

In this scenario the CRAB optimization works as fol- 
lows. For the guess control pulse d'Li(t) we take 

te [o,td] 

'Sxo)cos^ ^^^zjf-git) + Sxo, t e {td,T2] 

(11) 

where both g{t) and X{t) are only defined in the time 
interval {t^, T2], but the expression remains the one given 
in Eq. (j?]). Such a choice for d-i(t) ensures that it is 
always bounded by Sxo and xo- The control pulse rfi(t) is 
then the time inverse of d_i(t), that is, di{t) = d^i{T2 — 

As outlined previously, here the goal is not only 
the minimization of the overlap infidelity 1=1 — 
I (?/'g I "0 (r2 ) ) P , but also the minimization of the occupancy 
in the middle trap. To this aim, we use the following ob- 
jective functional, namely the cost function we want to 
minimize: 

J=l + ^ r dt[AEit) + \M^c,t)\^] 
^2 Jo 



T2 - 2td 



dt\{Mt)\m)f 



Ws\{iJsmT2))\' 



(12) 
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where |V'd) is the spatial dark state, which corresponds to 
the first excited state obtained by diagonahzing at each 
time the single particle Hamiltonian H{t) ~ /2m + 
V{x,t), AE{t) = \E2{t) + Eo{t) - 2Ei{t)\ with £;„(t) 
n = 0, 1, 2 the first three eigenvalues at time t of H{t) {Ei 
is the energy of the spatial dark state), and xc is the po- 
sition of the minimum of the middle well, where the node 
of the spatial dark state should be located. The energy 
difference AE{t) is used to keep the energy of the spatial 
dark state equidistant from the energies of the ground 
and second excited states, and reduce the transitions out 
of the dark state. The second line of Eq. corresponds 
to the projection of the evolved state on the actual spa- 
tial dark state in the interval [id,72 — id], whereas the 
weights We, Wg, and Wd can be adjusted for convergence. 
We note that a similar objective functional has been used 
in Ref. [H. 

We first investigate the behavior of the SAP process 
without optimization. In Fig. [7] we show the overlap fi- 
delity = 1 — I as a function of the transport time T2 
when the trap parameters v±i(t) are chosen to be time- 
dependent (blue-bright line), which fixes the minima of 
the triple well potential to the same energy level. Instead, 
the black-dark line corresponds to the scenario for which 
v±i{t) = lyt G [0,T2]. In this case the first three lowest 
eigenstates of the Hamiltonian H{t) are not in resonance 
and the evolved state tries to follow the second excited 
eigenstate, as also illustrated in Fig. [7] by the behavior 
of the overlap infidelity at long times. This phenomenon 
occurs because when the two outer wells approach the 
middle one, the minima of the outer wells correspond 
to a larger energy than the minimum of the middle well 
(see also the magenta dashed line in Fig. [2]), and there- 
fore it is energetically more favorable for the system to 
follow the second excited state. Additionally, we note 
that fixing the minima of the triple well configuration to 
the same energy level yields a spatial dark state whose 
node is not localized in the centre of the middle trap, as 
desired, but it is lifted towards the outer well with lower 
barrier (or potential maximum) . In order to have a node 
in the middle trap one should also require maxima at the 
same energy level, but this in not contemplated in our 
study as outlined in Sec. HT] Thus, the dark state has to 
be properly engineered and this also explains our choice 
for the cost functional (fT2|). 

As shown in Fig. [7] (bright-blue line on the top) , the 
asymmetry of the potential, due to the fact that the max- 
ima of the triple- well potential are not fixed to the same 
energy level (while the minima are), does not enable the 
atom to follow an eigenstate of H{t), in particular the 
spatial dark state, and an oscillatory behavior occurs. 
Interestingly, the occupancy in the middle trap, defined 
as 

Vc{t) = / dx\^{x,t)\\ (13) 

Pc{T2)^^ r dtr,(t), (14) 




0.2 0.4 0.6 0.8 1 




Figure 7. (Color online). Overlap fidelity = |(V'g|V'(r2))|^ 
(top) and probability of occupancy of the middle trap (bot- 
tom) vs. time. The blue (bright) lines refer to time-dependent 
trap parameters v±\{t), whereas the black (dark) lines refer 
to time-independent trap parameters v±\{t) = 1. The in- 
sets show zooms of both the overlap infidelity (top) and the 
occupancy probability (bottom). The diamond and square 
symbols on the top of the blue (bright) line of both panels 
represent the situations for which the CRAB optimization 
has been performed. In both cases Sxo = 1.43 /^m. 



is higher when we force the minima of the trapping poten- 
tial to be at the same energy level (i.e., time-dependent 
v±i). Here x™^(i) are the positions at time t of the 
maxima of the trapping potential of Fig. [1] (lower panel) . 

In figure [5] we show the optimal control pulses (i_i(t) 
[panel (a)] together with the corresponding probabilities 
of occupancy Vc{t) [panel (c)] for the transport time 
T2 = 31.4 ms, that is, the red diamond symbol in the 
insets of Fig. [7] (we recall that di is the time inverted 
control pulse of d_i). Instead, in panel (b) the trap 
parameter v-i(t) is displayed, which corresponds to the 
optimal quantum dynamics for T2 = 31.4 ms with time- 
dependent trap parameters and with objective functional 
given in Eq. (fT^ (wi is basically the time inverted of 
V-i). For the CRAB optimization we used Ng = 25 and 
we set Wg.d ~ 0.5 and we = 1- Such a choice is a good 
trade-off between small overlap infidelity and occupancy. 
In Table [H we show the results of both the overlap infi- 
delity and probability of occupancy, for which a further 
improvement was not possible. Indeed, our attempts at 
optimizing by considering d±i(t) as independent control 
pulses, by looking for an optimal set of frequencies Wfc, 
instead of coefficients {Ak,Bk}, or by varying td, have 
not been able to further improve the obtained results. 

Given our findings in Fig. [71 we have also analyzed the 
scenario for which v±i = 1. In this case we used for the 
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optimization the following objective functional 

J - wx [1 - |(V^g|V(T2))p] +ii;pPc(T2), (15) 

where the weights wj, wp are adjusted for convergence. 
The results of such analysis are illustrated in both TableU 
and Fig. m 

Concerning the final probability of occupancy in the 
middle trap, Pc{T2), the optimization decreases the value 
from 0.1166, obtained with the initial guess control pulses 
dj-i, to the value of 0.0466, obtained with the optimal 
control pulse, when the system tries to follow the second 
excited eigenstate; from 0.1498 to 0.0771 or to 0.0916 in 
the other two cases when the system is forced to follow 
the spatial dark state (see also TableU]). As illustrated in 
Fig. [8] (c), the probability distribution function Vc{t) is 
peaked around t = T2/2, that is, when the trap separa- 
tion is minimum. For comparison, we also show the prob- 
ability distribution for T2 = 44.8 ms (the green square 
symbol in the insets of Fig. [7|) , after optimization. In this 
situation, the distribution is almost symmetric with re- 
spect to t = T2/2, because the optimized dynamics tends 
to split the transport of the wave packet from the left 
to the right well not directly, but in a two-step-like pro- 
cess, where at time t = T2/2 the state is almost a dark 
state. This fact is also confirmed by the pair (pd,'T-d) 
of values of the projection onto the instantaneous spa- 
tial dark state [second line of Eq. (fT2|) ] and the node of 
the spatial dark state [second integrand in the first line 
of Eq. (HI])]: for T2 = 31.4 ms we have (0.211,0.027), 
whereas for T2 = 44.8 ms we get (0.151,0.011). Thus, for 
longer times, the system follows more closely the spatial 
dark state when the objective functional (|12p is chosen. 

A final remark on the choice of the objective func- 
tionals comes from the following fact. As shown in the 
second and third rows of Table U both the overlap in- 
fidelity and the probability of occupancy in the middle 
well are decreased when the objective functional ([TS]) is 
used in the CRAB optimization. This choice implies a 
better achievement of our goals as well as a computation- 
ally less demanding simulation. Indeed, with (|15p there 
is no need to diagonalize the instantaneous single particle 
Hamiltonian H{t). Hence, with this choice it is possible 
to efficiently optimize the transport in optical superlat- 
tices containing different interacting atomic species [g^l 
and, particularly relevant for our purposes, the transport 
of a condensate. In this case the determination of the in- 
stantaneous spatial dark state and its eigenvalue are more 
involved (65| . 

Finally, as in Sec. lIII A] we investigated the robustness 
of the optimized dynamics against trap position fluctu- 
ations by adding a shaking term to the control pulses 
d°^l{t), as in Eq. ([5]). Moreover, since the CRAB al- 
gorithm determines the optimal time-independent coef- 
ficients Ak and Bk, we were also able to investigate the 
effect of an imprecise control of the time delay td on the 
overlap infidelity. The result of such analysis is illus- 
trated in Fig. |9]for Wshake/w^: = 10~^ and for T2 = 31.4 
ms. We compare the obtained results for two cases: when 



72 (ms) 


I 




31.4* 


0.0007 


0.0466 


31.4'l^ 


0.0035 


0.0771 


31.4 


0.0048 


0.0916 


44.8 


0.0028 


0.0699 



Table I. Optimized errors for different transport times of the 
SAP technique for a single atom. The first row corresponds 
to the case for which v±i = 1, whereas in the other cases v±i 
are time-dependent. For the first two rows the objective func- 
tional used in the CRAB optimization is defined in Eq. (I15|) . 
but for the second row we fixed the minima of the trapping 
potential to the same energy level. Instead, for the last two 
rows, the objective functional is defined in Eq. (|12p . 



the atomic system is forced to follow the instantaneous 
spatial dark state (right) and when it is forced to fol- 
low the second excited state of the trapping potential 
(left). As already pointed out in Rcf. |38i], for purely har- 
monic confinement, the SAP process is less affected by 
imprecise timing, but the overlap infidelity drops faster 
for fluctuations in the trap positions than in the ideal 
case considered in Ref. [Sg. Indeed, while for harmonic 
traps the fidelity is reduced by 1-2% [sl] , at the optimal 
time delay and Oghake ~ ^x/2 {ix is the harmonic oscilla- 
tor ground state width), for Gaussian traps the fidelity 
< 60% (right picture), which shows how detrimental the 
anharmonicity of the trapping potential is [66j . How- 
ever, this phenomenon is also due to the fact that the 
system follows the second excited state, which is more 
sensitive to energy losses, and therefore to a worsening of 
the transport fidelity. Additionally, the figure shows the 
fragility of the dynamics that the instantaneous spatial 
dark state of the system is forced to follow, especially 
concerning the control of the time delay. This further 
analysis confirms what we already noticed in the overlap 
infidelity and occupancy probability, that is, it is more 
efficient to follow the second excited state rather than 
the spatial dark one. 



IV. TRANSPORT OF A CONDENSATE 

In this section we investigate the optimal transport of 
a BEC in optical dipole potentials such as the ones in 
Eq. ([1]). We assume the quasi-lD regime of quantum 
degeneracy and a mean field description of the atomic 
system dynamics, that is, we assume that the Gross- 
Pitaevsku equation (GPE) [67| 



d'ilj{x, t) 
dt " 



2m dx^ 



+ Vix,t)+giDN\ij{x,t)f 



(16) 
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Figure 8. (Color online), (a) Optimal control pulses obtained by means of the CRAB algorithm respectively for the three 
situations listed in the first three rows of TableU T^i =31.4* ms (blue dashed line), T2 =31.4 Vs (black solid fine), and T2 =31.4 
ms (red solid thick line), (b) Trap parameter v-\{C) for the optimal control pulse T-2. =31.4 ms (red solid thick line) of panel 
(a), (c) Probability of occupancy of the middle well vs. time obtained with the optimal control pulses for T2 =31.4*ms (blue 
thin line) and Ti =31.4 ms (red thick line), (d) Probability of occupancy of the middle well vs. time obtained with the optimal 
control pulse for Ti, =44.8 ms. Here the time delay is ~ 5.2 ms for Ti = 31.4 ms, td — 7.5 ms for Ti = 44.8 ms, and for both 
transport times 5a;o = 1.43 /im. 



well describes the physics of our problem. Here t/' 
is normalized to one, N is the atom number, g\D = 

2fe.^i.a|p (^1 - 1.4603^) [sll, is the three di- 
mensional (3D) s-wave scattering length, and aj^ = 
{2h/'muj±)^/'^ . This assumption implies that the radial 
confinement is frozen to its ground state, and therefore 
that the ratio 77 = uj^/lo^ is significantly larger than 1 
(a;_L is the transverse trap frequency, that is, the trap 
frequency in the y — z plane) . As we already pointed out 
in Ref. [2Q|, a good value is ry = 20, for which radial exci- 
tations due to two-body collisions can be suppressed. We 
underscore that a simulation of the current experimental 
setup would require a 3D simulation of the GPE, since 
the transport occurs in the transverse direction while the 
axial confinement has a shallower trap than the trans- 
verse one, where the transport process occurs. 



that for > A'cr the condensate is not stable and the 
GPE has no longer a stationary solution. We have in- 
vestigated this phenomenon in the quasi ID regime by 
considering the Gross-Pitaevskii energy functional. For 
an arbitrary confinement potential V{x) the functional is 
defined as |67| : 



E 
N 



dx 



2m 



dip{x) 



dx 



Vix)\^Pix)\' 



(17)" 



For a harmonic trap, a Gaussian Ansatz for the conden- 
sate wave function can be used, which has been proven 
to provide an excellent estimation of A^cr for a three- 
dimensional BEC [13 ■ To this aim, let us consider the 
following Ansatz for the condensate wave function (nor- 
malized to 1) 



A. Attractive inter-particle interaction 

It is well known (sec, for instance, Ref. [S^]), that for 
attractive interactions (i.e., gm < 0, that is, a^jj < 0), 
a critical number of condensed atoms exists, A^cr, such 



^(x) = {ae^V^r'/' exp (-^^j , (18) 

where cr is a dimensionless parameter which represents 
the effective width of the BEC. By inserting ([T51) into 
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Figure 9. (Color online). Transport efficiency (i.e., 1 — X) from (^^J(a;) to (fi^^\x) as a function of the time delay td between 
the two trap approaches and the amplitude Oshake of shaking the positions of the outer traps with tJshakc ~ W~^uJx- (left) for 
time-dependent trap parameters v±i, that is, when the system is forced to follow the dark state; (right) for time- independent 
trap parameters v±i{t) = 1, that is, when the system is forced to follow the second excited state. Here Sxq = 1.43 /im and the 
"optimal" time delay td — 5.2 ms. 



(fT7)) wc obtain 



E 



Nhw^ 



fcOt 



(TV27r 



(19) 



where /^^ = l-1.4603a|p/ai, and a = TjNal^/i^. The 
behavior of the GP energy functional, for some values of 
ij and N for an atomic cloud of ^^Rb atoms trapped in 
the hyperfine state |F = 2,mi^ = — 2), is illustrated in 
Fig. [TUl As the figure shows, the local minimum disap- 
pears for cither a small atom number (solid vs. dashed 
lines) or for a small ratio r/ (solid vs. dashdot lines). Im- 
portantly, the energy local minimum appears always for 
E < 0, that is, the interaction energy, Ei^t, exceeds the 
kinetic energy. Contrarily to the repulsive case, for at- 
tractive interatomic forces the kinetic energy, i^kim can- 
not be neglected. Indeed, it stabilizes the condensate 
against collapses, namely the condensate is stable as long 
as E'kin > Eint- We can give a rough estimation of A^cr 
by using this inequality: 



which imply 



2m£2 



E, 



int 



(20) 



(21) 



By fixing ry = 20, we have Na ~ 0.3, for utx = 2tt x 711 
Hz, and N^- ~ 2.7, for = 27r x 7.11 Hz. The latter 
value of Ncr could be in principle further enhanced by 
reducing ry, even though the quasi ID condition would 
no longer be very well fulfilled. Hence, from this analy- 
sis, we see that, for an attractive BEG in the quasi ID 
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Figure 10. (Color online). Gross-Pitaevskii energy functional 
vs. the effective width in the Gaussian model for attractive 
interacting atoms in the quasi ID regime. The black (solid) 
line corresponds to 77 = 20 and A'' = 10, the red (dashed) line 
to 77 = 20 and N = 5, and the blue (dashdot) line to = 5 and 
N = 10. For all curves we considered *^Rb atoms confined in 
a harmonic trap with frequency Ux = 2n x 711 Hz. 



regime, the admissible condensate atom number is ex- 
tremely small. With such condensate atomic numbers 
the realization of attractive BEGs in optical microtraps 
is actually not feasible. Given this, the attractive case 
will be discarded in the subsequent transport analysis. 
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B. Optimization of step 1 of the transport process 



Firstly, wc computed the ground state of the left well 
by using the imaginary time propagation technique. We 
start by considering 3N atoms in such a way that we have 
precisely N atoms in each well. This approach is valid 
when the three wells are far apart. The imaginary time 
propagation for 37V atoms yields a wave function which 
is the sum of three inverted parabolas, for large N, or 
almost Gaussian functions, for small N. Each of these 
three spatially separated density profiles is localized in 
one of the three wells. Then, we selected the part of the 
wave function localized in the left well, properly normal- 
ized, as initial state. This new wave function describes 
precisely a BEG with N atoms [we have checked the sta- 
tionarity of the solution by propagating it in real time via 
GPE with N in the mean-field potential and as confine- 
ment the initial (static) potential]. When the two outer 
traps move towards the middle one, since the Gaussian 
potentials do overlap, the three wells cannot be treated as 
independent anymore (see also Fig. [T] lower panel). Even 
though the condensate wave functions of each of the three 
wells do not overlap, the curvature of the outer wells is 
different from the middle one. Hence, it is no longer 
straightforward to determine how many atoms are con- 
tained in each well. To overcome this problem, we propa- 
gate adiabatically the condensate wave function trapped 
in the well centered in x'q towards the one centered in 
xq. The trap position xq is chosen as the minimum sep- 
aration between the three wells such that N\'tj}{x,Ti)\'^ is 
localized only on the left well, that is, the atomic density 
in the middle and right wells can be effectively neglected. 
We underscore that the value of xq crucially depends on 
N, for fixed gn). Hence, the adiabatically evolved con- 
densate wave function is chosen as goal state of step 1 
of the transport process. We then analyzed how fast the 
initial state can be propagated towards the goal state, by 
simulating the dynamics for different transport times Ti , 
which have been chosen much smaller than the adiabatic 
transfer time. 

We have considered atomic ensembles with = 10 or 
200 ^^Rb atoms per well (the latter have been recently 
obtained in experiments |5l|)- For 10 atoms the initial 
potential configuration is the same as for the single par- 
ticle scenario (Fig. [1] top), with xq = — 3.0/im and a 
slightly smaller trap frequency — 2tt x 702 Hz with 
respect to '^^Rb. For N = 200, since the condensate 
wave function has a much larger width than the single 
particle one, the dipole potential has to be adjusted. To 
this aim, in order to keep the (initial) lattice periodic- 
ity fixed, that is, \xq\ = 6.5//m (i.e., with laser beam 
waist 1.3 A<m), the potential depth has been increased up 
to Vo = fcs X 25 ^K. Such a potential depth implies a 
single-well trap frequency = (4Vb/m(T^)^''^ ~ 27rl2 
kHz, very similar to the trap geometry of Ref. [5l| . and 
Xq = — 3.8/im as minimal (target) separation. 
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Figure 11. (Color onlme). Overlap infidelity vs. transport 
time: black (solid) line with the control pulse -Dho(i) defined 
in Eq. ((22]); the other two lines show the infidelity obtained 
with the CRAB optimized control pulse dopt(i) = Dho{t)gopt{t) 
for Ng = 8 (dashed line) and Ng = 16 (dashdot line). The 
upper panel corresponds to N = 10, whereas the lower one to 
N = 200. 

As initial guess for the control pulse we used 

DUt) ^\vj,t-^ ^<t <^ , (22) 

where Aa; = \x'q\ — \xq\, and Vm = 3Ax/2Ti is the maxi- 
mum trap velocity during the transport. Such a control 
pulse has been proven to be optimal for a ID condensate 
in a moving harmonic potential at the transport times 
Ti,„ = 3{2n+l)T:/u}x with n e N [36]. Thus, there exists 
a minimum transport time, Ti^q = 3tt/uJx, for which no 
excitation in the condensate is produced. 

In exactly the same way as for the single-atom trans- 
port, we investigated the (quantum) speed limit of step 
1 of the transport process, whose results are illustrated 
in Fig. [H for iV = 10 (top) and N = 200 (bottom) 
atoms with repulsive interaction. For 200 atoms the po- 
tential depth is about Vq — 43.5huJx, the chemical poten- 
tial ^ = 39.2hujx whereas in the Thomas-Fermi limit we 
have fiTF = AhhJx [1^. Thus, the system is well within 
this limit. Goncerning the optimization, the GRAB al- 
gorithm works precisely as we described in Sec. IIII A[ 
with the only difference that we have to substitute the 
Schrodinger equation with the GPE and define the con- 
trol pulse as d{t) = Dho{t)g{t), where g{t) is given by 
Eq. (U). Besides, the overlap infidelity is defined again 
through Eq. ([2), where ^^'^{x) = ^{x,Tad) is the state 
obtained adiabatically. Tad 3 ms, starting from the 
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ground state of the left well of the initial potential con- 
figuration with trap separation |a;g| = 6.5 ^m. The same 
procedure is used for 10 atoms. 

We see from Fig. [TT]that, while for TV = 10 the infi- 
delity decreases monotonically with respect to the trans- 
port time Ti, for N = 200 this is not the case, and it 
becomes a monotonic function only for Ng = 16. We 
attribute this behavior of the infidelity to a non perfect 
revival of Bogoliubov excitation modes present during the 
transport process, which have a larger impact for bigger 
condensates, because of the larger non-linear interaction. 
To further improve the results one could also optimize 
the dynamics of the Bogoliubov collective excitations, for 
instance, by solving the time-dependent Bogoliubov-de 
Gennes equations [70|. This approach, however, would 
allow to engineer the Bogoliubov modes, but at the ex- 
penses of a very demanding numerical optimization. 

Furthermore, as shown in Fig. [TTl the transport times 
for N = 200 are shorter than in the single-particle and 
small condensate cases. This is basically due to a shorter 
transport distance Ax — 2.7 /im (in the single-particle 
case Ax = 4 /zm) and to the trap frequency, which is 

16 times larger than in the single atom scenario. Fur- 
thermore, we see that the control pulse ((22)) is an excel- 
lent guess with satisfactory overlap infidelities up to 1 
ms for 200 atoms and up to 16 ms for 10 atoms, even 
for transport times Ti 7^ Ti^„. Notably, with respect to 
the single particle optimization, the addition of harmon- 
ics does not improve significantly the overlap infidelity 
for short transport times. This behavior may also be 
related to the initial guess for the coefficients {Ak,Bk}, 
for which wc always started by setting their initial val- 
ues to zero. Indeed, this may occur also for Ti = 0.7 ms 
{N = 200), where the control pulse with Ng = 16 yields 
a slightly worst overlap infidelity with respect to the one 
obtained with Ng ^ 8. The choice for the initial values 
of {Ak,Bk} might be not the right one, since the control 
landscape may have several minima: the larger the num- 
ber Ng is. the larger the control landscape. Thus, our 
initial choice likely produces an optimal control solution 
trapped in a local minimum that is not the same for a 
larger Ng. We also note that by performing the opti- 
mization on the frequencies cuk instead of the coefficients 
{j4fc, Bk} the improvement in the infidelity is very small. 

Regarding the quantum speed limit [U, it can be 
roughly fixed to 0.5 ms for = 200, which is larger than 
l/i^x — 0.08 ms and is slightly smaller than Ti^o — 0.13 
ms. We (numerically) defined the limit by considering 
the time for which the infidelity is approximately 10~^. 
This is a reasonable threshold to quantify the error on 
the distance between the state evolved until time Ti and 
the goal state. We note that the quantum speed limit 
is roughly determined by maxt^[Q Ti]{h/[Eg{t) — Ee{t)]}, 
where h is the Planck constant, Eg(t) is the instanta- 
neous ground state energy, and is the instantaneous 
energy of the first excited state. For shorter times, it is 
not physically possible to bring the system in the ground 
state of the trap without populating excited energy lev- 
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Figure 12. (Color online). Difference between the guess con- 
trol pulse Dho{t) and the optimal one of step 1 of the transport 
process for Ti = 0.5 ms and Ng — 16 (A^ = 200). In the inset 
the optimal control pulse is displayed. 

els, which, on the other hand, are needed, during the 
interval (0,Ti), to perform a fast transport of the atom. 
Instead, for = 10 atoms, the quantum speed limit can 
be roughly fixed to Ti = 16 ms, where the infidelity is 
about 10""^. Even though in the single particle scenario 
we had a slightly higher value of the trap frequency, be- 
cause of the use of ®^Rb atoms, it is interesting to note 
that already a small atomic cloud significantly alters the 
quantum speed limit, which, for a single atom, has been 
estimated around 3 ms. 

In Fig. [T2] the difference between the initial guess 
and the optimal control pulse for N = 200 atoms, 
Ti = 0.5 ms and Ng ~ 16 is depicted. This plot shows 
how small is the correction on the guess control pulse, 
even though it is quite important to decrease by more 
than an order of magnitude the overlap infidelity. For 
= 10 atoms, the initial guess pulse (thick line in 
Fig. fO)) is rather different with respect to the CRAB op- 
timized one. This larger distortion is due to the fact that 
since the potential depth Vq = 2.55hojx and (initial) trap 
separation are the same as in the single particle scenario, 
the potential wells are not deep enough to consider the 
control pulse of Eq. (f22|) . optimal for a harmonic trap, as 
a good transport pulse. Indeed, while the single-particle 
energy is about OAQhojx, the chemical potential for 10 
atoms is /i ~ 1.87hujx- 

Finally, we also investigated the robustness of the op- 
timal control pulse for A^ = 200 against ffuctuations 
of the outer trap positions like for the single atom dy- 
namics. For the transport time Ti = 0.5 ms the op- 
timal solution obtained with CRAB is rather robust: 
the overlap infidelity changes from 0.0012 to 0.0046 for 
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Figure 13. (Color online). Control pulses of the step 1 of 
the transport process for TV = 10 atoms: initial guess given 
by Eq. (|22[) thick (black) line, and optimal CRAB pulse with 
Ng — 16 thin (red) line. The transport time is Ti ~ 15.9 ms. 



flshakc ~ — O.lfim. This effect is due to the coopera- 
tive behavior of the atoms in the collective motion of the 
condensate. Instead, we did not investigate the effect of 
dimensionality, because, unlike in the single particle sce- 
nario, the nonlinear term appearing in the GPE is also 
affected by the augmented space geometry, and therefore 
the comparison would not be fair (apart from the issue 
of validity in the quasi-2D regime). 



C. Optimization of step 2: SAP process 

The optimization of SAP with interacting particles is 
more difficult with respect to the single atom scenario. 
Indeed, as also discussed in Ref. (65l |. in the spectrum 
of the nonlinear Gross-Pitaevskii Hamiltonian loops 
near the avoided crossing points and new eigenstates of 
Hgp emerge when enhancing the nonlinear interaction. 
As pointed out by Graefe et al. [6^1, within a three-mode 
model, SAP, in order to work in the nonlinear regime, has 
to fulfill the following two conditions: (i) gioNA > 0; 
(ii) gioN/lx < gc = A. Here A represents a detuning 
between the three wells, that is, the resonance condition 
needed for SAP. We note that the resonance condition 
in this case imposes that the onsite energies of the wells, 
tvjjk{t) + fJ-kit), are constant at all times, where uJk{t) is 
the local frequency of the fc-th well and iikit) the corre- 
sponding chemical potential at time t. The inequality (ii) 
shows that there exists an upper bound on the nonlinear 
interaction strength for the realization of SAP. The prob- 
lem we are studying, however, cannot be strictly treated 
within a three-mode approximation. Nevertheless the 
model will be used as a guideline when discussing the 
results of the optimization. 

As for the single particle study, we applied the CRAB 
algorithm in order to understand whether optimal con- 
trol can improve the performance of the SAP protocol. 



Both for N = 50 and 200, however, we noticed that for 
a fixed number of harmonics {Ng = 10) CRAB was not 
able to reduce the value of the overlap infidelity obtained 
with the initial guess control pulses (jlip . This (empiri- 
cal) observation holds both when we are optimizing the 
control pulse by searching for the optimal set of coef- 
ficients Ak, Bk, and when we seek the optimal set of 
frequencies uJk- Moreover, we numerically noticed that 
the convergence of the algorithm to the value of the over- 
lap infidelity obtained with the initial guess control pulse 
takes longer than in the single-atom case. Even though 
in these two cases the number of atoms is likely much 
larger than the one allowed for the realization of SAP 
in the nonlinear regime, we attribute the occurrence of 
such a phenomenon to a more elaborated control land- 
scape topology, that is, a control landscape with a large 
number of local minima due to the emergence of new 
eigenstates in the system. We did not further investigate 
this aspect, which deserves a deeper analysis in a sep- 
arated work, but we rather chose to further reduce the 
number of atoms to iV = 10. In this case CRAB was 
able to improve the performance of the protocol with 
respect to the initial guess control pulse. As already 
pointed out in the previous section, with respect to the 
single-atom scenario, here we used *^Rb atoms which im- 
ply a smaller trap frequency lj^ and a trap separation 
xq = —3.0 fim. Apart from these small changes, due 
to a different atomic species and a broader size of the 
atomic sample, the trap configuration is essentially the 
one of the single-atom case. Nevertheless, the optimiza- 
tion carried out for different transport times T2 could not 
go below ^20% of overlap infidelity and ~10% of pop- 
ulation in the middle trap. The result of such a study 
is illustrated in Fig. [TH The obtained results cannot be 
improved by further optimizing the frequencies ujk ■ This 
shows that even though optimal control can improve the 
performance of the protocol, there is however a physical 
limit due to the SAP resonance condition for which no 
further optimized dynamics can be achieved. Indeed, at 
t = T2/2 the separation between the wells is minimal and 
we can roughly estimate the detuning as A ~ 0.15 hiOx, 
whereas gioN/ix — 5.56 Hujx, which shows how condi- 
tion (ii) is not satisfied even with only = 10 atoms. 
To increase A one should further reduce Sxq , but then the 
three wells merge in a single one, or, alternatively, by re- 
ducing the atom number. In this case, however, the BEC 
would be very small and the GP description might be 
also questionable. Although with a different trap setup, 
the analysis carried out in Ref. [131 also shows that the 
overlap infidelity increases quite quickly with the num- 
ber of atoms and that even with only two *^Rb atoms the 
(non-optimized) performance of SAP is quickly harmed 
(^16% of infidelity). Besides, as Fig. [HI illustrates, the 
behavior is not monotonic, which is probably related to 
a non optimal dynamics of the Bogoliubov modes. 

Finally, concerning the population of the middle trap. 
Fig. [U shows that it is almost constant with a mini- 
mum of about 0.1. We note that, in comparison with the 
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Figure 14. (Color online). Overlap infidelity (top) and prob- 
ability of occupancy of the middle trap (bottom) vs. time 
for the optimization of the SAP protocol for A'^ = 10 inter- 
acting *^Rb atoms. The minimum allowed trap separation is 



Sxo = 1.42 fj,m and v±i{t) — 1. 



single-atom case, we did not further minimize the pop- 
ulation of the middle well, since the transport efficiency 
was already lower, and therefore we preferred to focus 
on the minimization of the overlap infidelity [i.e., we set 
wp = in Eq. (jisp j. Nevertheless, the CRAB optimiza- 
tion has been able to further reduce the population with 
respect to the one obtained with the initial guess control 
pulse. 



V. CONCLUSIONS 

In this paper we have numerically investigated the per- 
formance of the SAP protocol by means of optimal con- 
trol both at the single particle and at the many-body 
level. In our analysis we have considered trap param- 
eters, atomic species, and atom numbers that are used 
in current experiments [7l] |. The transport process has 
been split in three steps, because of the initial large trap 
separation. The first step brings the atom(s) localized in 
the left well closer to the middle well in such a way that 
tunneling between the three wells occurs, therefore en- 
abling the realization of the second step of the transport, 
that is, the SAP process. Afterwards, the third step of 
the transport process brings further away from the mid- 
dle well the atom(s) localized in the right well. We have 
seen that while we can easily achieve the quantum speed 
limit, both for the single particle and the condensate sce- 
nario, for the first and last steps of the dynamical trans- 
port process, the second one requires a higher degree of 
control already for small transport time reductions with 



respect to the "adiabatic" times. In the single atom case, 
we observe a smaller population in the middle trap when 
the system is forced to follow the second excited state of 
the trap (i.e., time- independent v±i) rather than follow- 
ing the actual dark state (i.e., time-dependent u±i). In 
the latter case, due to the different energy level of the 
maxima of the triple well configuration, the node of the 
dark state wave function is not localized within the mid- 
dle trap, but outside. This fact forced us to additionally 
engineer the shape of the dark state wave function ren- 
dering the control landscape more complicated. Thus, 
we had to make a trade-off between transfer efficiency 
and suppression of the middle trap population. In addi- 
tion, we observed that the engineering of the dark state 
reduces the robustness against trap and time delay fluc- 
tuations of the optimal control pulse. We note that, in 
order to further improve the transfer efficiency and re- 
duce the population of the middle trap, by engineering 
properly the potential, one could employ a programmable 
and computer controllable nematic liquid-crystal spatial 
light modulator, where the trap separation can be varied 
by changing the periodicity of the modulator [t^I ■ Alter- 
natively, optical superlattices can be used, which would 
allow to fix the three minima at the same energy level as 
well as the two maxima. 

The optimization of the SAP protocol for a condensate 
strongly relies on the atom number and onsite energy of 
the wells. We have investigated in some detail the perfor- 
mance of the protocol for TV = 10 atoms with repulsive in- 
teraction. The analysis showed that the CRAB algorithm 
is able to improve the transport efficiency with respect to 
the one obtained with the initial guess control pulse, but 
the maximum attainable efficiency, for a transport time 
not longer than 450 ms, is about 80% with a population 
in the middle well of about 10%. It is not clear whether 
longer times could yield a better efficiency, which would 
require a longer computational time, but if this would 
be the case, one has also to take into account the effects 
of decoherence. For instance, if we consider atom chip 
technology 3 , where the expected limits due to surface- 
induced decoherence of motional states are comparable 
to the ones of the hyperfine states, which have coherence 
times of about Is 73[, our analysis already shows that 
we are actually close to the limit of the SAP protocol. 
This ultimate limit, for a relatively small BEC, is due 
to the emergence of new eigenstates and crossing levels, 
as discussed in detail in Ref. [gI] , which break down the 
SAP protocol when the nonlinear interaction exceeds a 
critical value. 

In summary, from our investigations, it emerges that 
while at the single atom level SAP can be optimized 
below the 0.1% level, and possibly observed in current 
experiments, the application of an optimized SAP tech- 
nique to a condensate is rather limited, already even with 
small number of atoms. On the other hand, it would be 
interesting to investigate more precisely and more gener- 
ally the influence of the nonlinear interaction of BEC on 
the quantum speed limit of a certain dynamical process. 
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and this will be pursued in future work. APPENDIX 
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The determination at each time of w-i(t) and vi(t) 
is a rather complicated nonlinear minimization problem. 
In our simulations, however, we noticed that an excel- 
lent approximation to the values of v^i(t) and vi(t) is 
given by the following procedure: at the beginning the 
positions of the minima of the trapping potential ([1]) are 
determined by looking for the roots {xl,xc,xr} of the 
function 



V'{x,t)= ^[x-fc4(t)]exp^ ^^ ^""f" , (23) 

fc=-l w ) 



where v±i = 1. Then, we use the following formulae: 
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These solutions are obtained by solving the system of lin- 
ear equations: V{xL,t) = V{xc,t), V{xji,t) = V{xc,t). 

Finally, we also mention that our numerical simula- 
tions of both the Schrodinger and the Gross-Pitaevskii 



equation have been performed by means of the split op- 
erator tech niq ue together with the fast Fourier transform 
algorithm [63|. 
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